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SUMMARY 



This investigation was performed at Princeton University in the 
period from Fehrviary, 1951 to June, 1951* 

The original intention was an experimental investigation of the 
feasibility of \ising a powered model in harnessed circling flight to evalu- 
ate the longitudinal stability of a full-scale airplane. The airplane 
selected was the Cessna l40 two-place light airplane, since its stability 
derivatives and other parameters had been determined and were at hand. A 
one-tenth scale flying model was constructed. 

The method of attack was to have been as follows: 

A. Determination of stick-fixed neutral point, from steady- 
state flight tests* 

B. Experimental determination of the model's frequency response- 
d 

-j by obtaining a transient response of 0 to an elevator pulse iiput 
and converting this to the frequency response by the method outlined in 
Reference 11. 

C. Correlation of A and B with theory. 

No quantitative data were obtained for A and B due to mechanical 
difficulties and lack of time. Sxnfficient qualitative data were obtained, 
however, to show agreement between the model's transient response and that 
predicted by theory; to point up severeLL difficulties peciiliar to this type 
of testing; and to Indicate recommendations for further development of the 
method. 

It was found that the scale chosen gave too small a model from con- 
siderations of required mechanical installations and dynamic similarity. 
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Transient response tests of relatively short dxiration were found 
to be the only feasible type due to the presence of vake disturbances. 
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IHTRODUCTION 



With the Increasing number of new and radically different configure* 
tlona of airplanes being Introduced for sonic, transonic and convertaplane 
purposes It seems that scmie method of evaluating the lov-speed longitudinal 
stability characteristics In unaccelerated forward flight Is desirable prior 
to fixing a design. This could most easily be done In model form. A tech* 
nique such as this would allow the Inertia forces on the model to be Inter* 
preted to full scale results and would also be valuable In determining gust 
loadings. 

The Air Materiel Command at Wright Patterson AFB established a test 
facility which gave promise of yielding such data. 

An analytical study (Beference 3) was made by AMC. From this It 
was concluded that certain of the longitudinal stability characteristics 
could be evaluated, and that even where the dynamic stability of the plix^gold 
mode was greatly modified by the harnessing forces on the model, the model 
response could be related to full scale response by certain mathematical 
operations. Unfortunately the test facility was decommissioned before any 
data vere taken to substantiate the above conclusions, and at the time of 
commencement of this Investigation, had not been re-established. 

This facility made use of powered models flown In harnessed circling 
flight; a type of flight known as "control-line flying" which is used as a 
sport or hobby by certain of the general public. 

It was the purpose, therefore, of the present Investigation to 
provide a limited amount of information on which to evaluate haraessed cir- 
cling flight of models as a testliig method, and to establish, If possible, 
testing techniques which will complement the method. To this end a one -tenth 
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scale model of the Cessna 140 light tvo-place airplane constructed. Reasons 
dictating choice of this airplane vere that complete static flight tests and 
vlnd tunnel tests of the Cessxia lUO hare been previously made here by the 
Aeronautical Engineering Department of Princeton University and therefore 
complete data vere available. This investigation vas limited to cciqparlng 
the neutral point, as obtained from level flight testing of the model, vith 
the neutral point from flight test of the full-scale airplane; and comparing 
the frequency response of the model's pitch angle to that obtained from the 
equations of motion of the airplaxie. The experimental frequency response 
was to be obtained from the analysis of the transient response to a triangu- 
lar pxilse of elevator deflection. 

These limitations vere dictated by lack of facilities, funds and 

time. 



The data vere to be obtained from static stability flight tests of 
the model with three center-of -gravity positions and from transient response 
teste at one center-of -gravity position (2J.d^ mac). 

It is believed that this investigation is the only one of its kind 
to date. Although the final results may not be conclusive due to rather 
primitive equipment, it may indicate a method and provide a part of the in- 
formation necessary for a complete analysis of the method, which of course, 
requires numeroxxs tests on many more of the longitudinal stability aspects. 
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AB 

AC 



a„ = 



a., = 



(— ) 

' dl o, / w 



at = 



b - 
c - 
C.G.- 

^Ow “ 
Cp = 

Ct- 

Cd: 



Wixig aspect ratio 
Aerodynaailc center of the ving 

- Slope of the lift c\irre of 
the ving (tvo dimensional) 

(per degree) 

- Slope of the lift curve of the 
vlng (three dimensional) 

(per degree) 

- Sloi>e of the lift curve of 
the tail (three dimensional) 
(per degree) 

Wing span — feet 

Mean Aerodynamic Chord of the Vlng — feet 



Mg 

Uo. 



Center of Gravity 
Dw 

D + Dw 



- Wire drag coefficient 

- Drag Coefficient 



Effective drag ooefficient (Cq-C-^) 



Parasite drag coefficient 



- Thrust coefficient 



C. 



'm 



Induced drag coefficient 

- -k. 

- b- 

- ^ - Pitching Moment coefficient 

c^Sc 

of the airplane 



- Lift coefficient 
~ Lift coefficient of the tail 
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Cio*. 




dL«. 




C 



- 



Azigle between flying 



i£‘- 



C 






d.«. 



C 







c 



'^de 




c 






= ^ 
d b 



Change in drag coefficient 
with change in angle of 
attack. 

Change In effectire drag coefficient 
with change in angle of attack, 
vires and datum plane 
Change in lift coefficient 
with change in angle of 
attack — per radian 
Change in pitching moment 
coefficient with change in 
cmgle of attack — per radian 
Change in pitching moment 
coefficient with rate of change 
of angle of attack with respect 
to 

Change in pitching moment co- 
efficient with rate of change 
of pitch angle with respect 
to 

Change in pitching moment co- 
efficient with change in 
elevator angle. 



D - 

d - 



e 



Total drag of the airplane — lbs 

Differential operator 

0.869 “ Airplane efficiency factor 
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- Non-dimens loml line force 



Centrifugal force in level flight 

£ £L 
R pvs 

Acceleration due to gravity — 32.2 ft/sec2 






- Non-dimensional inertia parameter 

- Airplane moment of inertia about 
the Y eucis — slug ft^ 



1 orj 



V^T 



mac 

- 

L - 
h - 



m - 

W\pk _ 




X - 

>u > 



V - 



Mean aerodynamic chord -- per cent 
Radius of gyration about the Y axis — ft 
Total lift of the airplane — lbs 

Tall length (airplane CG to AC horizontal tedLl) -- feet 

Moment about the Y axis — lb -ft 

Mass of airplane — slugs 

Miles per hour 

Stick-fixed neutral point 

Normal acceleration (along the Z axis at the CG) -- "G" vmits 
Period of oscillation — seconds 

p V ^ - Dynamic pressure 

Radius of fli^t circle. 

Wing area — ft^ 

p 

Horizontal tall area -- ft 

Time to damp to l/2 astplitude -- sec 

Time — sec 

Speed ratio , incremental change in airspeed divided 

V 

by airspeed 

True airspeed -- ft /sec 
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Volumetric coefficient of the tail 



W - Grose weight (average) -- pounds 
GRBEaC SIMBOIS 

^ - Angle of attack of vii]g — degrees 

cx^ - Angle of attack of the horizontal tall >- degrees 

- Elevator deflection — degrees 
£ - Angle of dovnwash — degrees 

~ - Airplane tall coefficient 

8 > Angle of pitch (angle between X axis and horizontal -- degrees) 
A - Root of stability quartic 

/Vy\ 

- — 5 — - Relative airplane density 

p Sc 

factor (non-dimensional) 
jO - Mass density of air -- slugs/ft^ 

T :: • ^ - r r - Time parameter — seconds 

({) - Phase angle — degrees 
to - Angular frequency — radians /sec 
Co - Angular frequency, radians per 



A variable with subscript zero indicates the maximum value of that 
variable during an oscillation and may be either positive or Mgatlve. 



Vertical bars indicate the conventional absolute valxie of any variable. 
SIGH COHVEirriOH 



Reference Axes 



Normal stability axes (wind axes) are used. 

X - axis - longitudinal axis - positive forward 
Y- axis - lateral axis - positive edong right wing 
Z-axls - vertical exls - positive down 
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Linear Dlsplacemeiits 

A linear displacwnent aJLong a positive reference axis is considered 

positive. 

Angular Displacements 

An angular displacement which is clockwise when viewed from the 
origin looking along a positive reference axis is coitsidered positive. 

Velocities and Accelerations 

Velocities and acceleration, either linear or angular, are con- 
sidered positive in the same sense as the corresponding displacements. 

Control Surface Deflections 

Positive elevator angle is associated with a downward movement of 
the elevator trailing edge. 

Phase Angles 

A negative phase angle signifies a time lag of the response to 
the forcing function. 

A positive phase eingle signifies a time lead of the response to 
the forcing function. 

ASSUMPTIONS 

In the theoretical analysis and in the evaluation of the stability 
derivatives the following assumptions are made: 

1. Stability axes (wind axes) are xised. 

2. The airplane is initially in level flight. ('Qxis means for 

stability sixes coi 6=1, Q - % ) 

3. All disturbances are small, so that the equations of motion may- 



be considered linear 
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The lift caused by the horizontal tail is negligible. 

Small lateral disturbances vhlch may exist do not affect 
the aerodynamic forces or pitching moments in the plane of 
symmetry. 

6. Forces due to pitching reloclties are negligible. 

7* Change in pitching moment due to the rate of change of 
elerator angle is negligible. 

8. Only longitudinal motion is considered. The moment of Inertia 
is that about the Y axis. 

9. Criteria for dynamic similarity were based on equal C i. of 
model and airplane (References ^ and 10). 

10. Standard model flight conditions used in the solution of the 
equations of motion were: 

V = 58.7 
B - 37.5 Jj^ 

at datum level = 4 ft. above the floor. 
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DESCRIPTION OF THE MODEL AMD IHSTRIMEiraATIOH 

The model used in this investigation vas a one-tenth scale model 
of the Cessna l4o tvo-place light airplane. It vas constructed by the 
authors. Die construction consisted primarily of balsa-wood, paper-covered 
and doped; a constiruction similar to spin tunnel test models (Reference 9)> 
Original Intentions had been to make the model dynamically similar on the 
basis of equal lift coefficient (References ^ and 10). However, the choice 
of airplane and the scale used soon Indicated that this would be physically 
Irqpossible with available equipment, due to the low weight requirement and 
small size of the model. The choice of airplane vas fixed by considerations 
previously mentioned, and scale vas fixed by (1) bxuigetaxy reasons plus 
(2) the requirement of an Indoor flying site to reduce stray wind disturbances 
on the model. It vas decided to stay as close as possible to the weight 
limitation and measure the resulting moment of inertia of the model, using 
that and the actual weight in the non-dimensional equations of motion of the 
airplane of the harnessed model. This modification in no way affects the 
ability to bhe effectlveiieas of this type of flight as a test method. 

The model vas x>overed by a Bantam .19> cubic inch displacement, 
l/lO HP (® 9000 RPM), Internal combustion, air-cooled, model airplane engine 
with glow ignition (Figure I ). This Ignition saved vei^t by requiring no 
ignition components in the model; merely a l|-volt battery as a starting 
device. The engine was mounted inverted to preserve the contours of the 
upper cowling and it projected out of the lover cowling a shor*t disteince 
(Figure 2 ). A pressiare type fuel feed system of commercial make was used 
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to prevent variation of engine speed due to change of fuel head In flight 
(Figure I ) . 

In order to give acess to the mechanism^ engine and fuel sTstem 
the vlng; as veil as various hatches vere made removahle (Figure ^ ) . 

In order to clarify the operation of the elevator pulsing mechan- 
ism It Is necessary to describe the control of the model In flight. In 
normal flight the model flies in a circle around the operator on two wire 

lines which lead frcaa a bellcrank (Figure ^ ) in the fuselage, through a gvdde 
To a 

near the wing tip^handle in the hand of the oi>erator. The model is flown 

with some yaw to the outside of the circle which is produced by offsetting 

the rudder and by attachliig the bellcrank behind the C.G. The amoxint of 

this yaw is also controlled by adjusting the line guide on the wlttg. A push 

«lft.vatoY Korn 

rod from the bellcrank links with tKe^so that a motion of the operator tend- 
ing to increase tension in the top line from the handle and decrease tension 
in the bottom line is transmitted as an up deflection of the elevator. 

The requirements on the elevator for the pxiLse input are that it 
(1) deflects from the trim elevator position, (2) returns to this position, 
and (3) remains locked during the ensuing transient motion. To this end a 
third wire from the operator to the model was employed. When a level equili- 
brium flight condition was attained, the operator trsinsf erred tension from 
the two flying vires to the third wire. The point of attachment of this 
wire was slightly aft of the model's C.G. in order to maintain slight outward 
yaw and hence line tension. The first result of tension on the third line 
was the application of a breike on the elevator pxish rod (Figxire ). This 
locked the elevator in the equilibrixaa position previously established. As 
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the tension was increased to counteract the full centrifugal force of the 
model, the third line also triggered a rubber-powered pulsing mechanism 
(Figure ^ ). The cam in this mechanism displaced the aft portion of the 
push rod, working against a spring which returned the elevator to the original 
locked position. The cam was cut to give a pulse of triangular shape with 
time, having a duration of approximately .04 seconds. By adjusting the posi- 
tion of the cam follower on the push rod, the aniplitude of the pulse was 
varied. The cam contacting the follower closed an electrical circuit throu^ 
a third wire which fired a flash bulb. After the elevator had been piilsed, 
the mechanism was automatically re-cocked for another run when the tension 
on the third line was released and the elevator again controlled by the ope3> 
at or. 

Measurement of elevator angle was necessary in both phases of the 
testing. In the steady-state tests the angle required for eq\illibrlvim at 
each velocity had to be known, and in the dynamic tests the determination of 
the exact size of the pulse was necessary. This angle was to be measitred 
with a Giannini Mlcro-torqxie. This instrument is an electrical potentiometer 
employing a rotating contact of very low mechanical torqiie gradient. Angular 
positions are measured electrically by suitable calibi^tion. In this applica- 
tion the microtorque, acting as a variable resistance, was mounted in the 
fuselage and attached by a linkage to the elevator (Figure H ). Eecording 
elevator position reqiilred twp leads from the model to the center of the cir- 
cle. Enameled steel flying wires were used for this piupose. From the oper- 
ator's position the leads were continued to the recordLing apparatus on the 
ground (Figure b ) . The steady-state elevator position was recorded by 
calibrating a voltmeter to elevator position. A recording oscillograph was 
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employed to record the piilse. 

No in-flight meoetirements vere possible in practice for tvo reasons. 
First, the vibration Induced by the running engine vas sufficient to destroy 
the electrical continuity of the micro-torque. Second, the friction intro- 
duced into the elevator system by the micro-torque linkage made control of 
the model marginal. 

Weights adjustable fore and aft vere used to alter the model C.O. 
position in running the steady-state phase. Normal C.G. for the model cor- 
responded to noxmal C.G. for the full scale Cessna and vas at 27. mac. 

This position vas used for the dynamic phases. 

Dynamic pitch-angle data vas recorded by means of a motion picture 
camera tracking the model from, the center of the circle. This camera vas 
mounted on a heavy tripod equipped vlth leveling bolts and a revolving platen. 

A flash bulb in the Sx'd line circuit, placed so that its fleush vo\ild shov on 
the film, vas used to mark the Instant of application of the pulse. The 
camera equipped vith a cross hair reticule. With the camera field verti- 
cal and the camera adjusted to track in the horizontal azimuth plane, these 
cross hairs vere a reference for both pitch angle measurement and vertical 
displacement. The camera speed vas calibrated for a time reference during 
the transient response. 

The use of an armory vas obtained for a testing site. This reduced 
the amount of stray air currents as ccmi>ared vith an outdoor site, and ln5)roved 
the chances of obtaining reproducible results. 
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TEST PROCEDURE 



Th© C.O. location of the model was found "by suspending it from 
three positions having vide angular separation about the pitch axis. The 
C.G. was at the intersection of plxmib lines hung frcm the x^oints of bus* 
pension. 

Tlie model’s moment of Inertia about the pitch axis was determined 
by swinging it as a coo^oxmd pendulvon. The method is described in Eeference 
6 Sind calculations given in Appendix G. 

This C.G. location mac) smd moment of inertia (.012 slug ft?) 

were used in csilculating parameters for the equations of motion of the model. 

Due to restrictions imposed by conditions at the flying site (low 
celling lamps) and by the light weight of the model, it was necessary to fly 
it from a 37*5 foot radius Instead of the 50 foot radius originally contma- 
plated. Ihls resulted in a higher rate of yaw, but enabled the outward yaw 
angle of the model ee set by the line guide on the wing, to be reduced be- 
cause of the increase in centrifugal force with decrease in radius. The 
shorter redlus was objectionable from the viewpoint of slipstream disturbances 
as the wake had less time to die out before being flown throiigh again by the 
model. 

A. Static Phase 

No static-phase runs were zoade due to the above-mentioned mechanical 
difficulties. For the determination of the stick-fixed neutral point the 
model wets to have been flown in steady level flight at a given C.G. position 
and velocity. The velocity would have been measured by timing a number of 
laps. The trim elevator position would have been measvured with the calibrated 
voltmeter. A number of flights were to be made at this C.G. position with 
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varlovis speeds. The speed vovild have been, varied by a combination of engine 
adjiistment and propeller change. This procedure of flying at various speeds 
vas to be repeated for a number of C.G. locations. 

B. Dynamic Phase 

With the pulsing mechanism vound the model vas flovn level and 
the speed vas measured. The electrical leads vere connected to the opera- 
tor's handle from the recording apparatus, and the recording osclUograidi 
started. The camera vas also started, and tracked the model under manual 
control. The pulse vas triggered, the pitch angle data being recorded by 
the camera and the elevator pialse by the oscillograph. The record vas con- 
tinued until the stable transient died out. 

RESULTS 

No static phase data vere possible because the mlcro-toz^jue angle 
Indicator voxild not record under the vibrating conditions produced by the 
engine of the model. 

A total of fifteen flights vere made during idilch the model vas 
adjusted for controllable flight and qualitative data vere observed for the 
dynamic phase. These consisted of approximate measurement of speed and 
visual observation of the following ; characteristics of the phvigoid mode, 
response to 5# 10, and 25-degree elevator pulses, and the behavior of the 
vake made visible by exhaust smoke. 

All elevator-locked flights vere above datxm. Speeds observed 
were between 30-50 mph. The characteristics of the phugold vere erratic. 

Its period was observed to vary betveen 1.5-2. 5 seconds. In some Instances 
It vas seen to diverge sharply. In one Instance, during 4o seconds of 
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elevator-locked flight, eeveral phugold oscillations vere seen to cornmence 
and das^ out, under the Influence of vake dleturhancea or other effects. 

The 25 “degree elevator pulse was found necessary for a visually 
observable pulse. 

At the last pulse during one of the runs, the energy remaining In 
the mechanism was low enoiJi^ to appreciably Increase the duration of the pulse. 
The result was a greatly Increased transient leading to a violent phugold 
oscillation which almost temlnated the entire experimental Investigation. 



DISCUSSION 

Static Phase ; 

Various power settings are required during teste for elevator power 
and stick-fixed neutral point \islzig harnessed flight. This results In changes 
In elevator power and the necessity for evaluating thrust outputs. Thus data 
reduction Is complex, and since these chazacterlstlcs may be found more easily 
from wind t\nmel tests it was concluded that their calculation frcm harnessed 
flight tests Is Inadvisable. Another deterrent Is external dlst\n:bance due to 
flight in the slipstream of the prevloiis lap. 



Dynamic Phase : 

The initial Intentions of this phase vere to collect and analyze 

the flight data xislng the photographs of the time history of the pitch angle 

0 

to plot — against time. This plot was then to be graphically analyzed to 
6 

find the frequency iresponse as described in Reference 11. This e:q)erlmentally 
determined frequency response curve was to be compared with the frequency 
response Ibom the analysis of the transfer function of the eqmtions of notion 
for the conditions of harnessed flight and of free fll^t. Since the equations 
of motion developed In Appendix A are non-dimensional and Include the harnessing 
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forces as a separate tens, they apply equally well to a model or to an air- 
plane dynamically similar. Because the small size of the model prevented 
the attainment of a dynamically similar moment of inertia, the equations, as 
solved, refer not to the full scale Cessna 140 but to eua airplane having edl 
the Cessna longitudinal stability derivatives and characteristics (Beferences 
7 and 8) except for the mass eind pitching moment of inertia. These two char- 
acteristics, found by direct measurement of the model and calculated, are 
those of an airplane dynamically similar to the model built for this investi- 
gation. 

It was assumed from the information in the AMC report (Reference 5) 
that the stabilizing and daiiq>ing effect of the harnessing line forces on the 
phugoid oscillation would eliminate this mode from any consideration in 
the analysis of flight data and woxild prevent it from becoming a problem in 
the testing techniqxie. This did not prove to be the case. Although on 
flights with the elevators lander control of the operator the model appesured 
stable with the C.G, at 27.85^ mac; when the elevators were locked the model 
exhibited a very short period phugoid (1^ to seconds period) when flown 
at elevations of about 10' above datum. Furthermore this phugoid oscilla- 
tion was quite unstable, increasing to double amplitude in less than one per- 
iod on numerous occasions. It was sufficiently unstable to make pulsing 
the elevator at the datum level (4' above the floor) when flying at 40 miles 
per hour possible only at the risk of destroying the model. It was decided 
that to enable data to be taken saJTely a higher datum plane was necessary, 
l.e., the operator standing on an elevated platform while flying the model. 
See sketch on page 19 . Since the location permitting this was not availabl , 
no quantitative data were obtained at present writing. Instead, the anphasis 
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of the InTestlgatlon was shifted to determination of the causes of the erron> 
eous prediction and to analyzing methods vhlch would make the fulfillment of 
the alms of the original Investigation possible. 




5ketch showing »v\etKod of 

dcxtum plaAS 

In an attempt to find the causes for the discrepancy between the 
observed flights and the theory, the assun^tlons used In writing the equa- 
tions of motion (Appendix A) wore re-examined. It was confirmed that the 
angle of the lines leaving the model In the drag direction could be neglected. 
This angle was less than 2^ under the most extreme flight conditions. It 
was also found that the angle In the gravity direction due to the wel^t of 
the lines was also negligible. Hence the assumption that the lines extend 
In a straight line from the operator to the model Is valid even In so light 
a model. This Is due primarily to the short radius employed. The conditions 
that the model must be charged with 3A drag of the harnessing wires 

and that the mass of the wires and their Inertia effects upon non-steady 
motions can be neglected (due to the small mass of the wire) were also veri- 
fied. Calculations svq>portlng the preceding are shown In Appendix D. 
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Eovever, obseznratlons shaved that at the shoirt radius flown (which was 
necessitated bj Indoor floor space and ceiling obstructions) and at the scale 
speed, the assumption that the circling flight would produce no aerodynasilc 
effects on the longitudinal motion was not completely justified. Calcula- 
tions show that the speed ratio between the velocity at the outer wing tip 
and the inner wing tip is 1.175 for the model used in this investigation. 

For comparison, to show that the short radius used here was no more detrimental 
to this condition than the radius used at the Wright Field testing facility, 
a calculation was made based upon their B-29 model (Reference 5)* This re- 
stilted in a speed ratio of 1.22 between outer and inner wing tip. Therefore 
the model used in the present investigation was less affected by this varia- 
tion of velocity with spanwise position than the large one flown at Wright 
Field on I50' lines. This speed ratio effect would produce a very complicated 
solution to the equations of motion. Therefore in the solution of the eqmtions 
of motion it is omitted, with the reservation, however, that at least portions 
of the discrepancies between the theoretical solution for the frequency response 
and the experimental must be chargeable to this omission. An addltioml aeio- 
dynamic effect produced by the circling flight path is the interference effect 
of the model's own wake. Flight observations showed an intermittent disturb- 
ance which was most easily noticed by a sudden inward banking oscillation of 
the model. This was determined to occur when the model's outer wing passed 
through the slipstream remaining from a previous circuit. The banking oscilla- 
tion was prodiiced by the fact that in counterclockwise circling flight the out- 
board wing was being influenced by an iqjward component of both the slipstream 
and the Inboard wing tip vortex. Confirmation of this reasoning was assisted 
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by the fact that the exhaust smoke of the engine made the slipstream vake 
visible. This vake vas seen to leave the model approximately tangent to 
the flight circle with some doimvash. See sketch belov. 




Ske-tch v'VCNke effects 

Since, at scale speed the model x>cksses the same azimuth point on the flight 
circle, and the vake appeared to widen as fast or faster than the ciarvatxure 
of the flight circle, the model crosses this turbulent area uiHess it is 
flovn above the previous path by a reasonable distance. 

All of the foregoing has not indicated any concltisive reason for 
the unstable phugold vhlch vas observed. It vas noted, hovever, that in 
both the ana lytical solution made by Braun In Beference 5 and in the IxrLtial 
solution for the present investigation, the familiar assvmrption of gliding 
flight condition, namely * 0 , vas used. The solution obtained for this 
condition (Appendix C) using the equations developed for harnessed flight 
developed in Appendix A, gave as roots for the modes of motion: 
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A.= -I.L38, -5.35 

- *001 ± .2631 



This indicates a very lightly damped, stable phugoid, vith a period of 5*1 
seconds vhlch danqps to 1/2 amplitude In 30 cycles. This, of course, vas 
coorpletely contreiry to experimental observations. It is also of interest 
to note that the short period mode has degenerated into two real roots. 

This Indicates that the model vas overdamped, caused by too low a moment of 
Inertia. This moment of inertia discrepancy had been expected. (See Des- 
cription of Model). 



characteristic eqmtions, both with and without the harnessing force, vere 
re-solved using the following thrust estimate; Cj * Cj) for all flight con- 
ditions close enough to the equilibria coxulltlon such that the small angle 



the chsu'acterlstic equation (Appendix C) yielded the following roots: 



This solution Indicates a trend In the right direction since It gives an 
unstable phugoid of essentlsuLly similar period to the glide condition (^.06 
seconds), but which increases to double its amplitude in 5*32 cycles. This 
corroborates the observed instability of the phugoid although it does not 
show the rapid divergence which was observed in tests. It also shows that 
the power has only a very slight effect on the short period mode. 

To check the validity of using the previovus thrust assumption an 
additional solution was made of the characteristic equation (Appendix C) 
assuming that vas constant at the value of at scale speed. This makes 



In order to evaluate the effect of thrust on the eqmtions, the 




-1.705 

r .0055 ± .265 A 



X„= -5.884 
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• # 

Cjj « 0 and = Cj) . The changes produced in the coefficients of the 
characteristic eqviatlon between this and the prevloxis assumption were so 

that they were within the accuracy of the calculation method used to 
factor the characteristic eqxiatlon. It was therefore concluded that the 
previous assumption was sufficiently precise for the size of the factors and 
for the method used in extracting the roots. 

As a further correlation of the period of the phugold obtained 
from the eqmtlons with the period from observations the value of the period 
was corrected for flight at condltioM other than at dat\im level. This was 
accomplished by the approximate relationship uj = . from 

V ^ 'fZ 

Reference 1 . Using the standard flight conditions and simple trigonometry, . 
a force of -.119 lUs per foot of height from datum was calculated, (measuring 
directions and forces in airplane stability axes direction). This is con- 
verted into a C - correction factor by dividing by the weight and adding 
Xs 

one; • Cj^ (1 -V* .078 H ). Therefore denoting by subscript ( condi- 

tions other than at datum, the following approximate relationship can be 
expressed: 




P 




Cu 

Cu(i + .078K) 



> 

] 4 - .078Vn 



Fran the above it can be seen that flying the model at I3 feet above datum 
approximately halves the period. Therefore the observed results agree 
reasonably with the analysis Insofar as the period of the i»hugold is concerned. 
Also the quality of the stability of this mode is in agreement, namely they 
are both unstable. The only discrepancy is the rapidly divergent oscillation 
which was observed but which failed to materialize frcmi the theory. This is 
unsjq>lained. 
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The foregoing discussion has been primarily about the phugold 
mode since this mode Indicated the discrepancy In the Initial assxm^tlons 
upon vhlch the equations of motion vere solved. Nov that reasonable agree- 
ment has been established betveen the theory and the qmlltatlve obseirva- 
tions, the asauB 5 >tlone which produced this agreement, namely , 

vlU be used In all fiurther emalysls. The phugold mode Is In Itself of no 
Interest In this investigation except as a complicating factor, since it Is 
most markedly affected by the harnessing force of the lines (Appendix C). 
Examination of the shoi^ period roots, on the other hand, show practically 
no difference betveen those of harnessed and Tinhamessed flight. It Is this 
fact that appears to make this testing method desirable for evEj.\iation of 
dynamic stability above the frequency of the phtjgoid. 

It Is considered Insufficient to compare the dynamics of complex 

systems on the basis of the characteristic equations alone. Such an analysis 

only considers the natural modes of motion of the system when dlstvirbed, but 

not under the action of a forcing fimction. Since the normal dynamic testing 

methods Involve the application of a forcing function. It is more reasonable 

to compare harnessed and free flight on the basis of the frequency response. 

This frequency response can be obtained in the conventional manner by solving 

0 

the transformed equations of motion for the transfer function — els shown 

o 

In Appendix B. The transfer functions to a step Input were reduced to the 
frequency domain by the graphical logarithmic method of Reference 2, eind 
illustrated in Appendix E. The frequency response, as obtained from the 
analysis of the equations of motion for both harnessed Eind unharnessed flight 
are plotted on Figures 7 ®rid Q . These are the cxjrves with which the 
experimental frequency response was to have been compared. 
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Examination of the frequency response curves Flgiures 7 and B 
show that above a frequency of 1.0 radian per the frequency response 

for harnessed flight is identical in shape vith that of the free flight, but 
the amplitude of the former is less. The fact that this difference in anpll-» 
tude is a constant factor is quite evident from the logarithmic plot of 
Figure 10 . On this plot there is a constant 2 decibel gain difference. 
Converting this gain difference into amplitude it was determined that the 
amplitude of the harnessed model was 8Cff> of the correct free fli^t valvie at 
all frequencies above Ci = 1.0. This loss of amplitude is not unreasonable 
since the lines offer restraining forces which would tend to limit the ampli- 
tude of oscillation in response to a forcing function. 

Since the aforementioned loss in amplitude is tied ig;) in numerovis 
terms of the transfer function, it would be difficult to evaluate it directly 
fixaa the haimessing force. However, it would appear that the change in am- 
plitude is a function only of the mass, speed, and flying radius of the model, 
^therefore it should be possible to evaluate this percentage change for a var- 
iety of permutations of these items and tabiilate the result to be used in ad- 
justing e:i^erimentally detemined frequency response curves to free flight. 

In order to see if the transfer function of the harnessed flight 

Q 

_ gave a transient response which was similar, qualitatively, to that 
b 

observed when the model elevator was piilsed; and also to establish desirable 
frame speed for the recording camera, it was considered necessary to solve 
for the transient response. In order to simplify this calcxilation, a square 
pulse was assumed instead of the triangular one actually imposed on the 
model. It shovild be noted that a triangular pulse was chosen for the model 
to enable the sin^lifled method of Reference 11 to be used in converting 
from the experimental transient response back to the frequency response. 
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The mechanical details of detemiixing the transient response from 
the transfer function are cdvered in the literature. Kvonerous methods may 
be used. The one chosen for use here is the Heaviside es^pansion which is 
described in Beferences 1 and li-, namely, 

Given the transfer function 



and roots of G(X) are X, , , Aj ' 

peated roots. 



^ Elii) 

A„ none of which are zero or re- 



Then §(t) = ^ (Ail = ^ 



C,(0) 



^ - I 



clX 



This solution is the time response to a unit step input (Appendl^c 
E). By using the method of Eeference 4 the response to a pxilse of finite 
width can be found. This method considers the pulse as being made of a 
positive unit step followed at time At by a negative unit step. By super- 
posing these two inputs a unit pulse of width At is obtained. By the prin- 
ciple of the Convolution Theorem the output is then the result of siQ>erposisg 
the outputs to these step inputs. These outputs are Identical in shape but 
of opposite sign and separated in time by an interval At. . These operations 
were accooiplished on logarithmic coordinates and the results are plotted on 
Figures 1 1 and I 2. . These are time histories of the pitch angle Q per 
unit elevator input. The solutions are for a pulse of .04 seconds, .10 
seconds, and .14 seconds. They show that the pitch changes of the transient 
are of reasonable amplitude for these two pulse widths, provided that large 
elevator deflections are \ised. This corroborates the vlsiially observed ef- 
fects. They also show that the pulse must be kept shoit or the transient 
response dxie to the short period mode is negligible conpared to the phugoid. 
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If a camera record at 100 frames per second were made, the pictures at each 
l/lOOth of a second would give ample data frxm which to plot an experimental 
curve similar to Figxire 12.. Also, since the slower decaying of the two short 
period roots was reduced to 1^ of its initial value at .7 seconds, only the 
first 7/lOths of a second would be used in the graphical analysis (Reference 
U) to obtain the fjrequency response of frequencies associated with the short 
period mode. 
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COHCLOSIOHS . 

The iiee of harnessed models in circling flight for evaluating 
longitudinal stability of full-scale airplanes has been investigated quali- 
tatively. 

It vas determined that this method is Inadvisable for finding elev- 
ator power and stick-fixed neutral point because of the variations required 
in test jwwer settings and becaiise of continuous slipstream disturbance. 

The method appeared favorable for obtaining transient pitch-angle 
response to a pulse elevator input, employing a high-speed motion-picture 
camera to record transients. 

There is a definite lower limit on the size of models which can be 
used in this method, considering dynamic similarity, due to the critical 
weight requirement and the i3ecessa3ry mechanical installations . 

The use of an indoor flying site feat^Iring a raised platform for 
the operator is recommended. This would permit initiation of transients at 
the datum plane of (flying wires horizontal) and would prevent destruc- 

tion of the model in the event of ensuing unstable phugoid oscillations. 

Previous to recording data the model should be flown below datum. 

To preclude effects of weike disturbance, data would have to be recorded with- 
in one lap of datum-level flight. 

If further investigation of this technique is intended with the 
Cessna configuration, it is recommended that a model of at least twice the 
size, i.e., one-fifth scale, be used. 
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APPENDIX A 



Development of genereil equations of motion to include harnessing forces. 
This development is essentially similar to that of Reference and uses the 
same assumptions. 

Assumptions are as follows. 

(1) Hie lines extend in a straight line from the operator to the 
model. 

(2) The height at which the model flies to maintain the harness- 
ing lines in the horizontal plane is referred to as datum height and 
the model flying at this height is flying in the datum plane. In 
this plane the wings of the model are assumed level. 

(3) The bank angle of the model (a<>) because of the action of 
the harnessing wire guide on the Inboard wing tip eqxials the angle 
the lines make with the horizontal. 




(4) For small displacements of the model from datum, 



1, and F remains constant. 

(5) The equations of motion use stability (wind) axes. 

(6) Aerodynamic effects due to curvature of the flight path 



are neglected. 
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The equations of longitudinal motion are as follows: 




= - Co' 1/9Y*S -VNsmif •*- CxV^pV'S 


(l-a) 


^Pe, =-Vr»Vir “ ' Cl. + Wcosa" aw V *■ s » v\ A (k 

' R 


(1-b) 


^ ■y — + ~ Cm ^ \J '"_S ^ 


(1-c) 



These equations are then linearised In the usual fashion by considering 

each term: oc. (D, •• 0 = ^ AkN/ Jr ^ F ^oc A ir + 'bF V 

■feV ?)cK. bck4 

After lotting Cq - Cj> “Ct this results in 
mV t(C*op 3 V)AV + (CdJ'/£|?V'' 5 w')ii^oi 4-Wzi-8’ = o 



^V^V» - (CupVS - z A<^) AV - ( Cu^ Vj^pV’-s'l Ao^ 

-Wim YAY -V /waN/'*'A^ = 0 

R 


(2-b) 


an/ " Cm ^ P c a ot 





~ = 0 (2-c) 

In order to put the term in terms of the longitudinal variables, 
substitute A((> = ^ (vhere is the vertical distance from datxm), and dif- 
ferentiate vlth resi>ect to time. Neglecting higher-order terms this gives: 



mVv - Waiww ay - Cup\/S AV - C u ^ o £» a (k -v f. ^ -O 

' R diX 

/ dlir = N/itvsn \mVn«v«. r = y. att 


(2-b)i 



Eovever, since initial conditions specify level flight at datum level, 
- O , and = Ay , and from equation (l-a), -VSs’\v\a($ 

therefore equation ( 2 -b) 2 ^ becomes 

WnVs ^ cv - CupVSaV - + ^ ^ ^ 
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Dividing tlirough equation a and b by s^£> and equation c by 

^pV*‘5c = %Sc and making eubstltutlons as follows: 



XK. 



- ^ 
V 
« 

- )L 

SI 



T - 



^ = 



pSM 

/v>S- 

aSc. 



r P ^ 

^ - ft "pTs 

^ 



and using the operator d to indicate differentiation with respect to 
and considering aXL •<.'>' as being the displacements from initial conditions 
(l.e.,^*,^«’ ) gives the linearized non-dimensional fona of the equation; 
Those are ejq>ressed in operator form. 



( <i + Cd^ + 



¥r 

AX I-' ^ 



1. Cs>„9L ■¥ X C u y - o 



Cudix + XC,i_^do<. -(dt^+lCi 
Z ^ 2 . 






"V = o 



(3-a) 

(3-b) 



-V Ld*” “ - JlCva*3o< 

^ \ 

+ S' - 

Vk 

In the analysis following, Co^ andC*^^ will be considered small and 
' neglected. 
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APPENDIX B 

Development of characteristic equation and transfer function ^ . 

Solution to the equations of motion (Eqviations 3^ Appendix A) is made 
by assuming a solution of form 



Ju. -=• >x, e 



At/r 



cx> 









vith the ° for the stick-fixed condition with A* measured from 

equilibrium elevator position. 

This resialts in 

(Co*' + XW, -V i(Cnc< ^ ^ 0* = 

X * 

C,_Aw., ( Cp* -V Cu<( ^ X t-plc(, -^/\-v X+-^^0, =0 (4-b) 

^ C W». «cl r X “ X ^ 0, = O (^-c) 

“vT K J 

For these homogeneous algebraic equations to have other than a trivial 
solution the determinants of their coefficients must vanish. 

Expansion of this 3rd-order determinant gives an equation of the form 



bjs \** V t>„X*^ V V bt X" -vb.X-vVo* ■=. O 



(5) 



By defining the combination » C d> s m Cr'’ , 

K 

the coefficients of the characteristic equation above are expressed as 



bc-- 



bM = 



= 

b. = 
W. ^ 

b„ ^ 



CxT - Cu., 

-ci ^ + ^ CuCp* - 

XX ^ ^ 

^ — Co^Chi.^- ^ 

X H 

+ -p(o 



•c.*) ■* Co*C..C 

1 

^ A- MCo*) 



e 

K 
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Solution for the transfer fianction ^ to a step input is made by 
using the Laplace transform on the equations of motion (Equations If). 

Since all initial conditions are assumed zero for the variables used, and 
the transform of a step function is unity in the Laplace transform as de- 
fined in Reference 3» this results in equations 6 below. 

(Cd -V V -5- Cl. G - o (6-a) 

2. i- 

CuXju V 4- Cuo i^ A-v-Pjo<. ^ ^ ^ 6=0 (6-b) 

^ ^ ^ r ^ ~ XIg “ ~_LCw\^i>e (6“C) 

h ~ ^ ~vT ^ 

Bote: Althoxigh the left sides of the equations above appear similar to 

eqxiatlons k, in equations 6, and 9 are variables of while in 

equations 4, and 6, are constsuits. Solving equations 6 by determin- 

ants for ^ gives a solution of the fora: 
oe 

^ 0-3 V ■•r Ck.X -V- 0.0 (7) 

b ^ V bqA'*-t baX’+'b^^C' A-b,X -V 
Where the b coefficients of the denominator are those of the character- 
istic equation and 

' = I 

z. 

Gk, " ^ C^Cuii - Ci_C C)^ 1 
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APPENDIX C 

Since the characteristic equation is the denominator of the transfer 

function f all solutions vill he made for the transfer function only, 

o 

All solutions vill be based upon the following so-called standard flight 
conditions (all angles in radians): 

V • 58.7 ft /sec ( mph) 

= .232 from ^ 

Cp * ,0616 from wind txuinel polar (Reference 8) corrected by 

addition of wire drag (Appendix D) 



w = 


1.53 lb 


R = 


37.5 ft 


• 


25.6 


s • 


1.59 ft2 


r s 


CVJ 

• 


K z 


.0843 


b . 


3.28 ft 


^ • 


.119 


s 


.0476 slug ft^ 


AR s 


6.76 


P • 


.00238 




«.4o6 from Reference 8 


C s 


.492 ft 




4.62 


e « 


.869 from Reference 8 



r - 



The following various thrust conditions were considered using harnessed 
flight equations: 

Condition I - C*r - O constant 

•Condition II - C-r-Cp at cvllot and xk (Cd^= o| 

Condition III - C, r - C p at equilibrium and aa. , and is 

constant at that value ( c - C 



p 



VP." — J 
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These condltioits produce the foUovrlzig chemges in the components of 



(1^;. and : 


Condition I 


Condition II 


Condition III 




1.0 


1.28 


1.28 


Cx>’‘ 


.0616 


0 


0 




.116 


0 


.116 




.0277 


.0269 


.0269 


Z 








The f oUovlng terms depend upon 


tail efficiency . The 


Condition values 


vere obtained from Beferences J 


and 8 and corrected for K* 


• 


Cv«^ j 


-.751 


-.964 


-.964 


Cw^cle 


-.149 


-.190 


-.190 




-.0465 


-.059^T 


-.0594 



Solution of the transfer functions for harnessed flight gave 

Condition I: &(\\ - -8.91 (A^ 2.37 X~ ,232 X ->• .003) 

(X* ^ 6.99 8.249 A* .7559A’- ^ .^55 A .018) 

Condition II: ® - -11.43 2.31 A t- .1437) 

^ ' ( A'* + 7.577 /V /10.I65 A' -V .4125 /\ ^ .702) 

Condition III: r -11.43 (A^-*- 2.31 A .1417) 

4 (X" -V 7 . 57 TA*+- 10.152 A'- + .3823 A *•- .702) 

Factored forms of the preceding are; 

Condition I; t(X\^ -8.91 (A v .0262)(A v .084 ) 2.236) 

* (\4. .0327)(A-»- l.i38 )(xZ 533)(X'-^ .2627!=!) 

Condition II: ^(>S ^ -11.43 ( X .063) (X 2.24^) 

6^ (Xi- 1.705 )(X-*- 5 .^)(K '-- .012 X -w .2&^) 

Condition III: Same as Condition II 
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J 

Condition II vas established as the most proper condition and the trans- 



fer function vas re-solved for unharnessed flight by ccasputing the coefficients 
for -P- = O . 

Condition Hi- 0 f x\ - -H.L3 (A -*- »012)(A v 2,298) 

Free flight: 6 (/\ i- 1.70lv)( A 5.&75)( - .0018 A -*■ .H3^) 
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APPENDIX D 



Calculation of drag effects upon the harnessing wires. This solution 
assumes no outward yaw of the model. Axes are as shown. 



Ay 




T - line tension • F “ 'TlJf 

Cjj * 1.2 two-dimensional flow Cp of cylinder at low "based 
on frontal area. 

Svs/ • area per foot of wire normal to velocity 
s .0007 for .008" wire 
= .001 for . 012 " \flre 

^ • 1/2 where at is velocity of wire at station x 




dD = = Kx’^dx 



Consider an element of wire, length <Ax , loaded by the drag Kx'’dx . 
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Taking the summation of forces In the drag direction, assuming T remains 
constemt over the element, 

2: Fy = Ts vw 0 - \< x’- - T" sU ( e V <16 \ - o 



T” i\v\ © ~ 1^ x'’ — "T" s » vs 0 <_o Jk dl 0 — V c,oSt 0 s <vs 1_© - o 

dl K ci X. 



eoS 0=1 , coS ci 6 ■= * ) S vvs I 0 -= •=. 

dx dx dx dx*- 



therefore 



\< x" - T c\\ 



- O 

d x'- 



Integrating gives 



43 - - 'Ijs^ ^ c 

dx '2.T 



Integrating again gives 

^=-'^ ^Cx-vO 

\t--r 

Solving hy constants for end conditions; 

T> - O 

C - '5ir^ 

Therefore iT-"r 

0] 

®l=i k - 3©], 



* K-^O 



? Vw 

For tension shared by all vires equally (most conservative condition) 



Operattrr 



®x^i " 

Drag of vires; 



(-13.72)(.0023^^) 57.3 * 1-8*^° 
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Control handle 
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D 



- 1 C^Deix - I 



r «25^*001) ^ S .0619 lb for .012" wire 

Dt*-»au " 2-Dwj^ z Ovj (‘L')3:.r)vj 

3 • 

- (3) ■ .144 lb. 

Detenalmtion of proportion of wire drag adding to airplane drag: 
Centroid of wire drag: 



Xo 




r^y. 

clD.ix 
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K x’dx 
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buLt 


D = < 

3 


^ 1 s 




Therefore 




\<X^ 

= H 


n.r 


3 V 






K x^ 
3 


-L.S 


h"L. 



IDierefore the wire drag centroid is located 3/4 of the distance along 
the radius from the operator. 

D wires • 3/4 (.144) « .108# 

Converting to conventional coefficient: 

Co^ « .108 = .0166 



Total model drag coefficient: 



Co --C 






= .040 .005 + .0166 = .0616 



"T o't A. V d V" ^ < 



D = Ct> £^S =. .OG»w(L.S2.) - 
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Calcxilatlon of Lina Tension; 



Centrifugal Line Force 



Side Force ® ^ * of 
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T, = ^ 

• »0^76 . (58.7)2 - 4.38 lbs 

37.5 



3_-C^=Cj^^(3 *.008 x 5 = .OLo 
' ~ “ (.o4o)6,52 - .26 lbs 



Total Itine Tension *4.64 lbs 
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APPEHDK E 



Frequency response ced.culatlon of model in harnessed flight, 






Components of the L»m-A*\^ plot are given Figure S and are added to 
obtain the cranplete frequency response in Figure 10 . This result is re- 
plotted on cartesian coordinates in Flgxires 7 emd S . 



‘ K = Z.'il 
Z6 T- 31 - 1 . S d t 



^ K.obSl V I 









<S 



^ CO 



- ^sCL^^' GO 















— ( — . I W*^ 



\ ^ 



(i^ \ 

\ .a LSI / 



1 



PRINCETON UNIVERSITY 
AERONAUTICAL ENGINEERING LABORATORY 



PAGE ^3 
REPORT I8l 



APPEHDIX P 



Tranelent response by Heaviside method: 

^ 

& GlCo'l ^ = 1 Q ( AO 

Harnessed flight, Ct = Co : 

e -ll.h3CA"-» 2.31 A .U39I 

6 ~ ^ 7.577 10.165 A^ + .1025 A + .702] 

P(o^ > -11. 43 (.1439) - -2.376 
Q ( o ) ” .702 

A, = -1.075: 

PC A.') = 2.907 - 3.939 .1459 

= -.886 

Q'CA.'^z 4 A,^+ 3(7.577) A'^ ^ 2(10.165) A, .4125 
= -19.82 66.079 - 3^.66 .f .4125 

= II. ot 

«< = .686(U.43) . -.7844 

A,Q’(A,') -1.075(12.01) 

- -5.885: 

PC AC) = 34.62 - 13.592 1- .1459 
= 21.174 

<^CAO= -814.8^-786.95 - U9.96 .4125 

" -147.40 

_ -11.43(21.174) = .2790 

AuOi’(Av')' -5.8w(l47.40) 

- -2.376 - .7844e"^*°’^^'^V .2790e5‘®Q*^*^V 

o 



(Short period mode only) 
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long period mode of harnessed flight, = C © ; 

= ( 7 it 135.5°) 

& 

The coD^ete solution for this condition is plotted in Figure 1 1 . 
Short period transient, free flight, Ct =C.d: 

(t\ - -2.383 - .559e"^*'^°^*‘'^- .2826e"^‘®^^*^'T' 
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APPEHDIX G 

Sample calculation of Moment of Inertia 

NN - Model veight - 1.53 11>. 

i - Pendulum length - 1.23 ft. 

^ : Swinging gear weight - .0331 Ih. 

Jii - Swinging gear pendulum length - 1.053 ft. 

T= Measured period of oscillation = 1.335 sec. 

V = Model volume ^ .216 ft.^ 

^ - Air density = .00238 slugs/ft.^ 

Gravitational acceleration = 32.2 ft 

sec. 2 

^ Inertia of swinging gear 
- ; .000975 

where t ® measiured period of oscillation of the gear alone 
Svibstituting in the above equation for ; 

I = .0124 slug ft.^ 
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1.31 sec. 
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Figs. 1 & 2 




FIOUrtE 1. - DISASSHI'IBLED KODEL 




FIOURE 2. - COMPLETE MODEL 



I 



i 



I 



Figs. 4a & 4b 
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figures La, Lb, - TOP VIEWS BEFORE ASSEPBLY SHOWING n’TRR'"AL 
KECHANISPiS 
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Figs. 4c ^ 4d 




••’IGUrtS l^c. - VIP'/ FROI-' RIGHT SIDS BEFORE ASSEMBLY 




FIGURE AB, - SHC'/ING PRIRARY CONTROL WIRES A'll CC’TROL HA'T)LE 



Figs. 4e ^ 4f 
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DETAILS OF HE/\3URiHG 
APPARATUS 
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Figure 8 




Figure II 
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